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BACKGROUND 


The  dry  deck  shelter  (DDS)  is  a  three-compartment  unit,  which  can  be  mated  to  a 
submarine  for  operations.  These  three  compartments,  the  hanger,  transfer  trunk,  and  chamber,  are 
described  below. 

(1)  Hanger:  With  a  dry  volume  =  1400  ft’,  it  is  the  largest  of  the  three  compartments.  It  has  a 
large  hatch  at  the  aft  end,  which  allows  for  the  ingress  and  egress  of  the  divers  and  their 
equipment.  The  hanger  can  be  pressurized  independently  and  flooded  with  water  to  various 
depths. 

(2)  Transfer  Trunk:  Connected  at  the  aft  end  to  the  hanger,  the  trunk  is  mated  to  the  submarine. 
Transfer  of  men  and  equipment  between  the  submarine  and  the  DDS  is  through  a  hatch  in  the 
floor,  and  may  need  to  be  done  when  the  hanger  is  flooded  and  open  to  the  sea.  Since  the 
submarine  is  maintained  at  a  pressure  of  1  atmosphere  absolute  (ATA),  the  trunk  is  often 
cycled  between  the  dry  (1  ATA)  state  of  the  submarine  and  the  wet  (pressurized)  state  of  the 
hanger.  Its  pressurization  and  flood  capabilities  are  the  same  as  the  hanger. 

(3)  Chamber:  The  chamber  is  attached  to  the  forward  end  of  the  trunk.  It  is  outfitted  for  DCS 
treatment  and  can  be  independently  pressurized  to  any  depth  required  for  treatment. 

During  operations,  varying  numbers  of  divers  are  present  in  the  compartments,  each 
producing  CO2  and  consuming  O2.  It  is  especially  important  to  control  the  level  of  CO2  in  the 
trunk,  where  there  may  be  7  or  8  divers  breathing  in  a  50-100  ft3  air  space  (when  the  trunk  is 
flooded).  Since  the  level  of  CO2  is  not  currently  measured,  this  control  is  accomplished  by 
periodically  venting  the  compartment.  Air  from  the  submarine  air  banks  is  let  into  the 
compartment  at  the  same  time  as  air  is  being  bled  out,  back  into  the  submarine.  This  procedure  is 
apparently  effective,  but  has  two  drawbacks:  it  requires  a  large  amount  of  air  and  is  very  noisy. 
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The  former  makes  keeping  the  air  banks  charged  difficult  and  the  latter  hinders  communication 
within  the  DDS.  For  both  reasons,  it  is  desirable  to  minimize  the  amount  of  venting  required. 

We  will  derive  the  equations  that  describe  the  amount  of  CO2  within  a  compartment  and 
then  apply  those  equations  to  the  problem  of  periodic  ventilation  of  a  closed  compartment. 
Finally,  we  will  compare  the  predictions  of  the  equations  with  some  measurements  of  CO2  levels 
in  a  DDS  from  actual  Fleet  operations. 


DERIVATION 


We  will  calculate  the  partial  pressure  of  CO2  in  a  given  compartment  by  determining  the 
number  of  moles  of  CO2  in  the  compartment,  assuming  that  the  ideal  gas  law  holds.  The  DDS  is 
operated  to  depths  of  165  feet  of  seawater  (fsw),  allowing  for  a  maximum  gas  pressure  of  6  ATA. 
This  is  within  the  range  of  validity  of  the  ideal  gas  law.  If  we  assume  that  the  respiratory  quotient 
is  unity,  then  each  mole  of  O2  used  by  a  diver  leads  to  the  release  of  a  mole  of  CO2.  Thus,  the 
only  change  in  the  total  number  of  gas  molecules  in  the  compartment  is  due  to  the  flow  of  air  into 
or  out  of  the  compartment. 


dntot  =  F 
dt  m 


=  F 


where  ntot  is  the  total  number  of  moles  of  gas  in  the  compartment,  and  Fin  and  Fout  are  the  the 
number  of  moles  per  second  flowing  into  or  out  of  the  compartment,  respectively.  The  net  flow 
into  the  compartment  is  F. 

Carbon  dioxide  levels  in  a  compartment  change  due  to  three  different  mechanisms: 
introduction,  production,  or  elimination.  Carbon  dioxide  is  introduced  by  the  inflowing  air  from 
the  submarine  air  banks,  which  contain  air  with  an  upper  limit  of  1250  ppm  CO2  (the  submarine 
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is  usually  ventilated  using  surface  air  before  the  air  banks  are  pressurized,  so  the  actual  CO2 


concentration  is  often  substantially  below  this  level).  Produced  by  divers  that  are  present  in  the 
compartment,  CO2  is  eliminated  by  venting.  (We  do  not  consider  the  elimination  of  CO2  by 
scrubbing.  This  technique  is  used  on  occasion,  and  only  in  the  chamber  portion  of  the  DDS.) 
Thus,  we  can  write  for  the  number  of  moles  of  CO2 


dn 


cm 


dt 


C-  ■  F-  +  N  j- 
in  in  diver 


•A  — 


n 


cm.F 

ntot  out 


where  nco2  is  the  number  of  moles  of  C02  in  the  compartment,  Qn  is  the  mole  fraction  of  C02  in 
the  inflowing  gas,  Ndiver  is  the  number  of  divers  in  the  compartment,  A  is  the  average  rate  of 
production  of  CO2  by  a  diver  (in  moles  per  second),  and  the  ratio  nco2/nt0t  mo^e  fracti°n  of 
C02  in  the  gas  in  the  compartment.  This  final  term  implies  that  the  gas  in  the  compartment  is 
well-mixed. 

There  are  four  cases  of  interest  (Some  sub-categories  may  also  be  explored  below). 

a)  Fout  =  0,  F  =  0  =>  Fin  =  0 

(constant  pressure  in  the  compartment  with  no  venting) 

b)  Fout  —  0,  F  0  =>  F^  ^  0 

(pressure  in  the  compartment  is  increasing) 

C)  Fout  >0,  F  =  0  =>  Fin  =  Fout 

(venting  the  compartment  at  constant  pressure) 

d)  F0ut  F  ^  0  =>  Fjn  ^  F0ut 

(venting  the  compartment  while  changing  the  pressure) 

In  order  to  simplify  the  discussion  below,  we  make  the  following  definitions: 


xi  =  ntot 
X2  =  nco2 

B  =  Cin*Fin  +  Ndiver*  A 
C  =  FoUt 


The  above  equations  then  become 
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dx^  1 
~dt 


=  F 


dx~  x-> 


Case  a) 


Fom  =  Fin  =  F  =  0  =>  C  =  0 


0=>  Jtj  =  Xjq 


dx^ 
dt 

dXry 

-3f=B^x2  =  x20  +  B-t 

where  Xio  and  X20  are  the  initial  values  of  Xj  and  X2,  respectively.  Thus,  the  total  number  of  moles 
of  gas  is  constant  and  the  number  of  moles  of  CO2  increases  linearly  in  time.  Since  Fjn  =  0,  the 
rate  of  increase  of  CO2  depends  only  upon  the  number  of  divers  and  their  rate  of  CO2  production 

(B  =  Ndiver*A). 


Case  b)  Fout  =  0  =  C,  F  =  Fjn  ^  0 

dx-. 

— ~  —  F-  =>  X,  =  X*  +  F*  ’ t 
dt  in  **1  ^10  in 

(lXr}  f  \ 

~dt  =  B=>  x2=  x2q  +  B-t  =  x2q  +  \C^n  ■  Fm  +  Ndiver  ■ Aj-t 

Thus,  with  only  inflowing  gas  both  the  total  number  of  moles  of  gas  and  the  number  of  moles  of 
CO2  increase  linearly  as  a  function  of  time. 

Casec)  Fout  >0,  F  =  0  =>  Fin  =  Fout 
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J 


dx- 


1  _ 


dt 

dx 2 
dt 


=  0: 


=  B 


X i  *10 


...2- . q  x  -B_.x 
Xir\  ^  W  A2  £  a10 


10 


1-e 


*10 


+  x. 


20 


This  is  the  case  where  the  compartment  is  being  ventilated.  Thus  the  mole  fraction  of  CO2  in  the 
compartment  exponentially  approaches  the  ratio  B/C,  where 


,j  C--F-+N • A 
B  _  in  in  dvver 

C~  F 

rout 

The  higher  the  rates  of  production  and  introduction  relative  to  elimination,  the  higher  the 
asymptotic  level  of  CO2 


Case  d) 


Fout  >0,  F  ^  0  =>  Fin  ^  Fout 


dx, 

-±  =  F=$xl=x10+F-t 

dXn.  Xn.  C'Xr, 

_ Z  D  Z  D _ Z 

dt  c  \ 


We  can  solve  this  equation  for  X2  by  changing  variables.  First,  let  t 1  =  Xio  +  F*t,  so  that  dt 1  = 
F*dt,  and 


B’-C- 


x2 

f 


where  B 1  =  B/F  and  C '  =  C/F.  This  equation  is  now  separable.  If  we  do  a  change  of  variable  a 
second  time,  defining  v  =  X2/t 1 ,  then  this  equation  becomes 


df  _  dv 
f  ~  B-(1+C>v 
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Consider  two  cases,  Fin  =  0  and  F;n  >  0. 

Subcase  dl)  F;n  =  0  =>  C '  =  C/F  =  F0Ut/(Fin  -  Fout)  =  -1,  so  that 


df_  _  t dv_ 
f  ~  B ’ 

or 

v  =  5’ln(0  +  v 


0 


Substituting  for  v  and  t '  from  above  yields 


xn 

no~ 


X,  r&F-t  ~  F  '  ln(x10  +  F  • +  V0 


(lln^o  +  FO  +  Voj-Uio  +  Ft) 


At  time  t  =  0,  X2  =  X20,  so  that 


*2  = 


F  t 
\  out 


ln( 


no 


- )  +  -20 

*10  Fout't  X10J 


(*10  F out  *) 


Subcase  d2)  Fjn  >0,  F  ^  0,  so  that 


df  _  dv 

f  ~  B-a+cyv 

MO  =  ■  ln(S-(l  +  C)  •  v)  +  D0 


where  Do  is  a  constant  of  integration.  Substituting  in  v  =  X2/t '  and  rearranging  gives 


.f-Dn.F-C 


2  1 +  C  0 


At  time  t  =  0,  t '  =  Xjo  ,  X2  =  X20  >  and 


1  =(-ML.x  -x  ).x  c ' 
0  \1+C”  “To  20  /  To 


Then 


B ’ 


'4 


B ’ 


2  "  1+C  *  ll+C’  'TO  A20 


C 
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Upon  substitution  back  to  the  original  variables 


out 


X2=~F-i-xV)-(Fin-Fout),)->--t- *10  - ' *20> 


in 


in 


10 


,  xin  -(F.  -F  )-t 
v  10  in  out  ) 


F.  -F 
in  out 


where  B  =  Ndiver*A  +  Cm*Fm  .  Note  that  this  case  corresponds  to  the  simultaneous  venting  of  the 
chamber  while  the  pressure  is  changing  (assuming  that  the  volume  of  the  compartment  is  fixed). 
The  partial  pressure  of  CO2  can  now  be  determined  using  the  ideal  gas  law 


R  T 

P  =— - x 

C02  V  2 

where  Rg  is  the  gas  constant,  T  is  the  absolute  temperature,  and  V  is  the  gas  volume. 

VENTING 

During  operations,  venting  usually  takes  place  when  other  factors,  Ndiver ,  V,  and  P  (total 
gas  pressure),  are  being  held  constant.  Thus,  between  vents,  when  there  is  no  flow  into  or  out  of 
the  compartment,  Case  a)  above  describes  the  buildup  of  CO2  in  the  compartment.  During 
venting,  the  inflow  and  outflow  are  balanced  (constant  P),  and  Case  c)  applies.  Ventilation  occurs 
periodically,  repeating  after  time  tcycie ,  with  ventilation  lasting  for  time  tvent-  Define 


E  = 


1  -E 


■t 


vent 


Then  the  number  of  moles  of  CO2  at  the  end  of  a  vent  is  given  by 


nfe  Nj.  ■ A  +  C .  F. 

oft  _  diver _ m  in  r^Yon  r 

X2  ~  F  10  *  2  '* 

out 
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where  the  superscript  indicates  that  the  venting  gas  was  turned  off  or  on  at  that  point  in  time. 
During  the  nonventing  period,  Case  a)  applies.  Therefore,  at  the  end  of  a  non-venting  period  we 
can  write 


x°n  =  xf  +N  • A-(t  .  -t  ) 

2  2  diver  cycle  vent 

We  can  combine  these  equations  to  generate  recursive  relations  describing  the 

CO2  level  at  the  beginning  and  end  of  each  vent  cycle  in  terms  of  the  levels  of  the  previous  cycle. 

So, 


~  N,.  ■ A  +  C .  F. 

off  dryer _ 

2,i  "  F  , 

out 

N,.  -A  +  C.  F. 
on  diver _ in  in 

2,i 


ijL — i-IL-xm-E,+\  x°ff  ,+N,.  ■  A-(t  j  -t  )]■£ 

10  l  2,1  — 1  diver  cycle  vent  J 


^,n-£’+4"  ,-E  +  Nd,„r-A-(  t 


out 


TO 


2,i  - 1 


cycle  vent 


) 


In  order  to  simplify  this,  let 


N..  A+C.  F. 
G  =  E’*XiU*  diver _ in  in 

F  + 
out 


and 


Then 


H  —  Ndiver*A.*(tcycle  _  tvent)* 


G  +  E  H  +  E-x°ff  , 

2,i  —  1 


G  +  H  +  E  ■  <  •  ! 

2,i  - 1 


Both  of  these  equations  are  of  the  form 


fi  =«•/,-!  +b 
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Which  can  be  rewritten  in  terms  of  the  initial  value  fo 


We  can  determine  the  asymptotic  value  for  this  expression  by  letting  the  number  of  cycles  go  to 


infinity.  Because  E  <  1,  a1  =  E1  -40,  and 


off  _G±EJL 

2,oo  ~  i-E 

on  G  +  H 
2-  i  -E 


the  difference  between  these  two  is  simply  H  =  Ndiver*A*(tcycie  -  tvent),  which  is  the  amount  of 
CO2  produced  by  the  divers  during  the  nonventing  period. 

Using  this  expression  for  X2,»on ,  and  assuming  we  know  the  upper  limit  of  CO2  partial 
pressure  (either  0.015  ATA,  or  possibly  0.02  AT  A),  we  can  determine  the  venting  time  necessary 
to  maintain  CO2  partial  pressure  below  this  value  by  solving  for  tvent  in  the  transcendental 
equation 


(1  -E)-R 


diver 


A  +  C. 
in 


diver 


cycle 


-t  J 

vent 


R  T 

_g _ 

V 


where  P2,~on  is  the  CO2  partial  pressure  limit. 

Appendix  A  contains  the  source  code  for  program  TABLE.FOR,  which  uses  the  above 
equation  to  solve  for  venting  times  in  each  of  the  4  gas  volumes  present  during  DDS  operations. 
These  4  gas  volumes  are  the  chamber,  the  trunk,  the  hanger,  and  the  bubble.  The  bubble  is  the 
forward  portion  of  the  hanger  that  remains  gas-filled  when  the  rest  of  the  hanger  has  been  fully 
flooded.  It  contains  the  controls  for  venting  the  hanger,  and  is  separated  from  the  rest  of  the 
hanger  by  a  plexiglass  partition  that  is  attached  to  the  upper  portion  of  the  hanger.  Also  included 
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in  Appendix  A  are  a  sample  input  file  for  the  program,  and  4  sample  venting  tables  generated  by 
the  program.  Note  that  the  program  uses  an  IMSL  FORTRAN  subroutine  DZBREN  to  find  the 
venting  times  from  the  transcendental  equation.  Other  equation-solving  routines  can  be  used  in 
its  place. 

Appendix  B  contains  the  source  code  for  program  MODEL.FOR,  which  uses  the  entire 
set  of  equations  above  to  model  expected  CO2  concentrations  during  the  operation  of  a  DDS. 
Also  included  are  sample  input  and  output  files,  and  some  comparisons  to  data  from  actual  DDS 
operations.  Note  that  this  program  has  been  designed  to  be  used  only  for  the  transfer  trunk,  since 
this  compartment  is  the  one  for  which  the  most  detailed  operational  data  is  available. 
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APPENDIX  A.  Source  code  for  program  TABLE.FOR 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cThis  program  will  generate  DDS  venting  tables  for  the  four 

cdifferent  areas  of  the  DDS;  the  hanger,  the  bubble,  the  trunk, 

cand  the  chamber.  It  is  assumed  that  the  air  in  each  of  these 

ccompartments  is  well  mixed.  It  is  further  assumed  that  the  trunk 

cis  half  flooded,  gas  space  volume  =  120  cubic  feet,  the  hanger 

cis  unflooded  during  venting,  gas  space  volume  =  1400  cu  ft,  the 

cbubble  is  half  flooded,  gas  space  volume  =  120  cu  ft,  and  that 

cthe  chamber  is  unflooded,  gas  space  volume  =  170  cu  ft.  The 

ctemperature  is  290  K  =  63  F,  and  the  flow  rates  at  different 

ccompartmental  pressures  are  38  SCFM  from  0-9  fsw,  160  SCFM  between 

clO-29  fsw,  and  240  SCFM  at  greater  depths.  The  input  parameters  are: 

c 

c  ichamb  =  compartment  of  interest 

c  1  =  hanger  output  in  HANGER.TAB 

c  2  =  trunk  output  in  TRUNK. TAB 

c  3  =  chamber  output  in  CHAMBER.TAB 

c  4  =  bubble  output  in  BUBBLE.TAB 

c  A  =  oxygen  consumption  rate  in  liters/minute 
c  pco2max  =  C02  limiting  pressure  in  ATM 

c  cppm  =  C02  concentration  in  ppm  in  venting  gas 

c  tcycle  =  venting  cycle  time 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

implicit  real*8  (a-h,o-z) 
real*8  k,ndart,nda 

dimension  vtable(4,20),depth(4),imin(4,20),isec(4,20),flows(4) 
external  dzbren,func 
common/comf^cinf,dinf,k, tcycle 
c 

c  set  constants 
c 

pi  =  4.d0*datan(l.d0) 

ftm  =  2.54d0*12.d0/100.d0 

ftm3  =  ftm*ftm*ftm 

rgas  =  8.31d0  !  units  of  Joule/mol*K 

temp  =  290.d0  !  in  K 

p0=l.d0  !  atm 

rgt  =  rgas*temp 

pOpa  =  1.01325d5  !  one  atm  in  units  of  pascal 
c 

c  get  parameters  for  run 
c 

open(unit=3,  file=fable.inp’,  type=bld) 
read(3 , 1 00)ichamb,a,pco2max,cppm,tcycle 
close(3) 

100  format(il2/8(fl2.7/)) 

co21imit=pco2max 
c 

c  convert  to  proper  units 
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c 

pco2max  =  pco2max*p0pa 
aaa=a 

tcycle  =  tcycle*60.d0 
c  convert  a  from  liter/min  to  mol/sec 

a  =  a/(rgas*temp)/60.d0/1000.d0*p0pa 
c 

c  set  op  arrays  for  printing  of  table,  depths  are  from  current  tables 
c 

depth(l)=9.d0 

flows(l)=38.dO 

depth(2)=29.d0 

flows(2)=160.d0 

depth(3)=100.d0 

if(ichamb.eq.l)depth(3)=66.d0 

flows(3)=240.d0 

depth(4)=165.d0 

if((ichamb.eq.l).or.(ichamb.eq.4))depth(4)=130.d0 

flows(4)=240.d0 

c 

c  convert  compartment  volume  to  cubic  meters,  set  maximum  number 
c  of  divers  for  the  compartment,  and  open  the  output  file 
c 

if(ichamb.eq.l)  then 
vol  =  1400.d0*ftm3 
ndmax  =  20 

open(unit=4,  file=hanger.tab’,type=’new) 
else  if(ichamb.eq.2)then 
vol  =  120.d0*ftm3 
ndmax  =  8 

open(unit=4,  file=frunk.tab’,type=’new) 
else  if(ichamb.eq.3)then 
vol  =  170.d0*ftm3 
ndmax  =  8 

open(unit=4 ,  file=  ’chamber . tab  ’, type=  ’new  ) 
else  if(ichamb.eq.4)then 
vol  =  120.d0*ftm3 
ndmax  =  4 

open(unit=4,  file=’bubble.tab’,type=’new) 
else 

write(6,110) 

110  format(’  improper  compartment  number  -  stop  program) 
stop 
end  if 
c 

cppm  =  cppm/l.d6 
errabs  =  O.dO 
errrel  =  l.d-6 
do  10  id=l,4 
depthdds=depth(id) 

c  convert  flows  from  cu  ft/min  to  mol/sec 

flow=flows(id)*ftm3*p0pa/rgt/60.d0 

c 

c  calculate  ambient  pressure  (in  pascal),  total  moles  of  gas  in 
c  the  compartment,  and  the  time  constant  for  the  compartment 
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c 

pamb  =  (depthdds/33.dO  +  I.d0)*p0pa 
xlO  =  pamb*vol/rgt 
k  =  flow/xlO 
c 

do  20  nd=l,ndmax 
ndart  =  nd*a*rgt 
nda  =  nd*a 
c 

c  calculate  coefficients  for  equation  in  the  DDS  report 
c 

cO  =  (1  -  E)  cinf  +  dinf*(tcycle  -  tvent) 
c 

c  where  E  =  dexp(-k*tvent) 
c 

dinf  =  ndart/pOpa/vol 

cinf  =  (pamb*(nda+cppm*flow)/flow  -  pco2max)/p0pa 
tlow  =  O.dO 
thigh  =  l.d4 
maxfn  =  10000 
c 

c  dzbren  is  an  EMSL  FORTRAN  root  finding  routine 
c 

call  dzbren(func,errabs,errrel,tlow, thigh, maxfn) 
vtable(id,nd)=thigh 
imin(id,nd)=thigh/60.d0 
isec(id,nd)=thigh-60.d0*imin(id,nd)+l 
20  continue 

10  continue 

volf=vol/ftm3 
jmax=ndmax 
if(ichamb.eq.  l)then 
write(4, 1 8 1 )  volf  ,aaa,co21imit 

181  format^  HANGER,  unflooded,  volume  =  ’,f6.1,’cu.  ft.  7 

1  ’  02  consumption  =  ’,f5.2,’liter/min7 

2  ’  C02  limit  =  ’,f6.3,’  atmT) 

else  if(ichamb.eq.4)then 
write(4, 1 84)volf,aaa,co21imit 

1 84  format(/’  BUBBLE,  water  depth  =  4.5  ft,  volume  =  ’,f6. 1 , 

1  ’cu.ft.7 

2  ’  02  consumption  =  ’,f5. 2,’ liter/ min’/ 

3  ’  C02  limit  =  ’,f6.3,’ atmT) 

else  if(ichamb.eq.2)then 
write(4, 1 82)  volf,aaa,co21imit 

1 82  format(/’  TRUNK,  water  depth  =  4.5  ft,  volume  =  ’,f6. 1 , 

1  ’cu.  ft.’/ 

2  ’  02  consumption  =  ’,f5.2,’liter/min7 

3  ’  C02  limit  =’,f6.3,’ atm’/) 

else  if(ichamb.eq.3)then 
write(4, 1 83)volf,aaa,co21imit 

183  format(/’ CHAMBER,  unflooded,  volume  =  ’,f6.1,’cu.  ft. 7 

1  ’  02  consumption  = ’,f5.2,’liter/min7 

2  ’  C02  limit  =’,f6.3,’ atm’/) 

end  if 

if(ichamb.eq.  1 )  write(4, 191) 
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if(ichamb.eq.2)write(4,190) 
if(ichamb.eq.3)write(4,190) 
if(ichamb.eq  .4)  write(4, 1 94) 

190  formats#  DiversU6x,DEPTH(FSW)7 

1’  0-9  10-29  30-  100  101  -  1657, 

2  ’  38SCFM  160SCFM  240SCFM  240  SCFM7) 

191  format(/’#  DiversM6x#DEPTH(FSW)7 

1’  0-9  10-29  30-  66  67  -  1307, 

2  ’  38SCFM  160SCFM  240  SCFM  240  SCFM7) 

194  format(/’#  Divers ’,16x,DEPTH(FSW)  7 

1’  0-9  10-29  30-  100  101  -  1307, 

2  ’  38  SCFM  160  SCFM  240  SCFM  240  SCFM 7) 

write(4,200)(j,(imin(l,j),isec(lJj),l=l,4),j=l,jmax) 

200  format(20(2x,i3,5x,4(i3,2h,  ,i2,lh”,5x)/)) 

stop 
end 
c 

real*8  function  func(x) 

implicit  real*8  (a-h,o-z) 

common/comf/c,d,xk,xc 

func  =  (LdO  -  dexp(-xk*x))*c  +  d*(xc  -  x) 

return 

end 


End  of  souce  code  for  program  TABLE.FOR 
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Sample  input  file  TABLE.INP  for  program  TABLE.FOR 


3 

1.25d0 

0.015d0 

1250.d0 

5.d0 

read(3 , 1 00)ichamb,a,pco2max,cppm,tcycle 
End  of  sample  input  file  TABLE.INP 
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Table  1.  Sample  Output  File  HANGER.TAB  From  Program  TABLE.FOR 


Depth  (FSW) 

0-9 

10-29 

30-66 

67-130 

Diver  No. 

38  SCFM 

160  SCFM 

240  SCFM 

240  SCFM 

1 

0’34" 

0’  13" 

O’ 15" 

0’  32" 

2 

r  7" 

O’ 25" 

O’ 30" 

r  3" 

3 

1’41" 

0’38" 

O’ 45" 

1’35" 

4 

2’ 15" 

0’51" 

1’  1" 

2’  7" 

5 

T  48" 

r  4" 

r  16" 

2’ 38" 

6 

3’ 21" 

ri6" 

1’31" 

3’ 10" 

7 

3’ 54" 

1’29" 

1’47" 

3’ 40" 

8 

4’ 27" 

1’42" 

2’  2" 

4’ 11" 

9 

4’ 58" 

1’55" 

T 18" 

4’ 40" 

10 

5’ 29" 

2’  8" 

T  33" 

5’  9" 

11 

5’ 60" 

2’21" 

2’  48" 

5’ 36" 

12 

6’ 29" 

2’ 34" 

3’  3" 

6’  3" 

13 

6’  58" 

2’ 46" 

3’ 18" 

6’ 28" 

14 

T  25" 

2’ 59" 

3’ 33" 

6’ 53" 

15 

7’ 52" 

3’ 12" 

3’ 47" 

7’  16" 

16 

8’  18" 

3’ 24" 

4’  x" 

7’  39" 

17 

8’ 43" 

3’ 36" 

4’  15" 

T  60" 

18 

9, 

3’48" 

4’ 29" 

8’ 20" 

19 

9’ 30" 

3’ 60" 

4’ 42" 

8’ 40" 

20 

9’ 52" 

4’  12" 

4’ 56" 

8’ 58" 

Hanger,  water  depth  =  0.0  volume  =  1400.0 

O2  consumption  =  1.25  liter /min 
SCFM  =  standard  cubic  feet  per  minute 
Tcycle  =  5  min 
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Table  2.  Sample  output  file  TRUNK.TAB  from  program  TABLE.FOR 

Depth  (FSW) 


Trunk,  water  depth  =  0.0  volume  =  120.0;  02  consumption  =  1.25  liter/min  Tcycle  =  5  min 

Table  3.  Sample  output  file  CHAMBER.TAB  from  program  TABLE.FOR 


Depth  (FSW) 

0-9 

10-29 

30-100 

101  - 165 

Diver  No. 

38 SCFM 

160  SCFM 

240  SCFM 

240  SCFM 

1 

O’ 35" 

O’ 13" 

O’ 24" 

O’ 48" 

2 

1’12" 

0’28" 

0’51" 

1’41" 

3 

1’50" 

O’ 43" 

1’20" 

2’ 34" 

4 

2’ 28" 

O’ 60" 

1’50" 

3’ 23" 

5 

3’  5" 

1’  18" 

2’ 20" 

4’  4" 

6 

3’ 38" 

1’37" 

2’ 48" 

4’ 38" 

7 

4’  8" 

1’56" 

3’  14" 

5’  6" 

8 

4’ 35" 

2’  15" 

3’ 37" 

5’ 28" 

CHAMBER,  unflooded,  volume  =  170,0  cu,  ft.;  O2  consumption  =  1.25  liter/min;  CO2  limit  =  0.015  atm  Tcycle  - 


5  min 

Table  4.  Sample  output  file  BUBBLE.TAB  from  program  TABLE.FOR 


Depth 

(FSW) 

Diver  No. 

0-9 

10-29 

30-100 

101  - 130 

38 SCFM 

160  SCFM 

240  SCFM 

240  SCFM 

1 

O’ 36" 

O’  14" 

O’ 25" 

0’35" 

2 

r  15" 

O’  29" 

O’ 54" 

ri5" 

3 

1’55" 

O’ 47" 

1’26" 

1’  58” 

4 

2’ 36" 

r  7" 

2’  1" 

2’41" 

BUBBLE,  water  depth  =  4.5  ft,  volume  =  120.0  cu.  ft.;  O2  consumption  =  1.25  liter/min;  CO 2  limit  =  0.015  atm 


Tcycle  =  5  min 
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APPENDIX  B.  Source  code  for  program  MODEL.FOR 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

This  program  will  model  the  C02  concentration  in  a  dry  deck 
shelter  (DDS)  under  the  following  assumptions: 


c 
c 
c 

c  1 

c 
c 

c  2 

c 
c 

c  3 

c 
c 
c 

c  4 

c 
c 
c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


parameters  include  the  pressure  in  the  compartment, 
the  depth  of  water  in  the  compartment,  and  the 
number  of  divers  in  the  compartment, 
fixed  parameters  include  the  02  consumption  rate, 
the  C02  level  in  the  venting  gas,  the  initial  C02 
level  at  time  t  =  0. 

the  maximum  flow  rate  is  depth  dependent 
38  SCFM  from  0-9  fsw 
160  SCFM  from  10-29  fsw 
240  SCFM  >30  fsw 

changes  in  compartment  volume  or  pressure, 
or  venting  of  the  compartment  do  not  occur 
simultaneously. 


implicit  real*8  (a-h,o-z) 
real*8  nda,nco2,ngas 

cccc 

dimension  tmin(500),fsw(500),dwater(500),ndiver(500),ivent(500) 
dimension  vol(500),pamb(500),dpambdt(500),tsec(500),dvdt(500) 
dimension  fin(500),fout(500),nco2(500),ngas(500),pc02(500) 

common/ddsO/rtrunk,pi 

common/comvent/pambO  1 0,pamb  1 03 0,  ventO  1 0,  vent 1 03 0,  vent30 

xnil=l.d-10 

zed=0.d0 

open  (unit=  1  ,status=bld’,name=  ’model. inp) 
open  (unit=2,status=’new,,name=’model.ouf) 

rg=82.0562d0  !  gas  constant,  (cc-atm)/(gmol-degK) 

rgas=8.31  !  gas  constant,  J/(gmol-degK) 

rgl=rg/1000.d0 
rgm=rg/l.d6 

rgf=2.89787d-3  !  gas  constant,  (ft3-atm)/(gmol-degK) 

temp=290.d0 
tempO=273.16dO 
rgast=rgas*temp 
ftm  =  2.54d0*12.d0/100.d0 
pi=4.d0*datan(  1  .d0) 
c 

c  estimate  the  radius  of  the  trunk  by  using  a  volume  of  240  cu  ft 
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c 

rtrunk  =  (240.d0*3.d0/4.d0/pi)**(l.d0/3.d0)*ftm 

pOpa  =  1.01325d5  !  latm  in  pascals 

p0=  l.dO  !  1  atm  pressure  for  reference 

fswO=33.dO  !  latm  in  fsw 

t=0.d0  !  initializing 

pambO  10=1 O.dO/fs  wO*pOpa 
pamb  1 030=30 .dO*pOpa 
ventO  10=38  .d0/60  .dO/rgf/tempO 
ventl030=160.d0/60.d0/rgf/temp0 
vent30=240.d0/60.d0/rgf7temp0 
c 

c  read  in  parameters 
c 

read  (1,*)  nnode  !  number  of  nodes  in  profile 

read  (1,*)  a  !  oxygen  consumption  rate  (liter/min) 

read  (1,*)  cin  !  concentration  co2  in  ppm  in  venting  gas 

read  (1,*)  cout  !  initial  cone.  co2  in  ppm  in  chamber 

a  =  a/(rgl*temp0)/60.d0  !  convert  to  mol/sec 

cin=cin/l.d6 

cout=cout/l.d6 

c  fvent=fvent/60.d0/rgf/temp0  !  change  from  SCFM  to  mol/sec 
c 

c  will  assume  that  venting  lasts  for  the  entire  period  from  the  current 
c  node  to  the  next  node,  assuming  that  venting  is  off  by  the  next  node, 
c  this  requires  that  when  ivent(k)  =  1,  tmin(k+l)=tmin(k)+tvent 
c 

do  10  i=l, nnode 

read(  1,11  )tmin(i),fsw(i),dwater(i),ndiver(i),ivent(i) 

pamb(i)=(fsw(i)+fsw0)/fsw0*p0pa 

dwater(i)=d  water(i)  *ftm 

vol(i)=vtrunk(dwater(i)) 

tsec(i)=tmin(i)*60.d0 

10  continue 

11  format(3fl0.2,2il0,fl0.2) 
do  12  i=  1, nnode- 1 

fin(i)=zed 
fout(i)=zed 
delt=tsec(i+ 1  )-tsec(i) 
c 

c  assume  that  chamber  volume  and  pressure  do  not  change  simultaneously 
c  nor  do  such  changes  occur  during  venting 
c 

c  calculate  fin  and  fout  for  changes  in  ambient  pressure 
c 

dpambdt(i)=(pamb(i+ 1  )-pamb(i))/delt 
if(dpambdt(i).gt.  1  .d-8)then 
fin(i)= vol(i)  *dpambdt(i)/rgast 
else  if(dpambdt(i).lt.(- 1  .d-8))then 
fout(i)=  -vol(i)*dpambdt(i)/rgast 
end  if 
c 

c  calculate  fin  and  fout  for  changes  in  volume 
c 
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d  vdt(i)=(  vol(i+ 1 )- vol(i))/delt 
if(dvdt(i).gt.  1  .d-8)then 
fin(i)=pamb(i)*dvdt(i)/rgast 
else  if(d  vdt(i)  .It.  (- 1  .d-8))then 
fout(i)=  -pamb(i)*dvdt(i)/rgast 
end  if 
c 

c  calculate  fin  and  fout  for  venting 
c 

if(ivent(i).eq.  l)then 
fin(i)=fventfcn(pamb(i)) 
fout(i)=fin(i) 
end  if 
12  continue 
tout=0.d0 
delt=60.d0 
c 

c  calculate  the  co2  partial  pressures  at  each  of  the  nodes 
c  in  the  profile 
c 

ngas(  1  )=pamb(  1 )  * vol(  1  )/rgast 
nco2(  1  )=cout*ngas(  1 ) 
do  200  leg=l,nnode-l 
nda=ndiver(leg)*a 
fint=fin(leg) 
foutt=fout(leg) 
fhet=fint-foutt 
cnet=foutt 
bnet=cin*fint+nda 
delt=tsec(leg+ 1  )-tsec(leg) 
xlO=ngas(leg) 
x20=nco2(leg) 
c 

c  cases  (a)  &  (b):  fout=c=0, 
c 

if(dabs(cnet).lt.  1  .d-8)then 
ngas(Ieg+ 1  )=f  1  case(x  1 0,fnet,ddt) 
nco2(leg+ 1  )=f2casea(x20,bnet,delt) 
end  if 
c 

c  case  (c):  fout=c>0,  f=0  ->  fin=fout 
c 

if((dabs(cnet).ge.l.d-8).and.(dabs(fnet).lt.l.d-8))then 

ngas(leg+l)=xl0 

nco2(leg+l)=f2casec(x20,xl0,bnet?cnet,ddt) 
end  if 
c 

c  case  (d):  fout=c>0,  fin  >  0,  f  <>  0 
c 

if((dabs(cnet).ge.l.d-8).and.(dabs(fnet).ge.l.d-8))then 
if(dabs(fint).lt.  1  .d-8)then 

ngas(leg+ 1  )=x  1 0-cnet*delt 
nco2(leg+ 1  )=f2cased  1  (x20,x  1 0  ,bnet,cnet,delt) 
else  if(dabs(fint).ge.l.d-l)then 
ngas(leg+ 1  )=x  1 0+fnet*de!t 
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nco2(leg+l)=f2cased3(x20,xl0,bnet,cnet,fnet,delt) 

end  if 
end  if 

200  continue 

do  201  i=l,nnode 
pco2=nco2(i)*rgast/vol(i) 
tout=tsec(i)/60.d0 
pgas=ngas(i)*rgast/vol(i) 
c 

c  convert  to  atm  and  ppm  for  output 
c 

pco2=pco2/p0pa 
pgas=pgas/pOpa 
cppm=nco2(i)/ngas(i)  *  1 .  d6 
c 

c  want  cppm  as  surface  equivalent 
c 

cppmsev=nco2(i)/ngas(i)  *  1  .d6*pgas 
write(2,l  1  l)tout,pco2,pgas,cppmsev 
write(6,l  1  l)tout,pco2,pgas,cppmsev 

201  continue 

1 1 1  format(3fl  1.5,3x,fl  1.2410) 

30  stop 

end 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccc 

real  *8  function  vtrunk(d) 
implicit  real*8  (a-h,o-z) 
common/dds0/rtrunk,pi 

v  trunk  =  piy3.d0*(2.d0*rtrunk-d)**2.d0*(rtrunk+d) 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccc 

real*8  function  flcase(xlO?flow,delt) 

implicit  real*8  (a-h,o-z) 

f  1  case=x  1 0+flo  w*delt 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccc 

real  *8  function  f2casea(x20,slope,delt) 

implicit  real*8  (a-h,o-z) 

f2casea=x20+slope*delt 

return 

end 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccc 

real*8  function  f2casec(x20,xl0,b,c,delt) 

implicit  real*8  (a-h,o-z) 

expterm=dexp(-c/xlO*delt) 

f2casec=b/c*x  1 0*(  Ld0-expterm)+x20*expterm 

return 

end 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CCCCCCCCCC 

real*8  function  f2casedl(x20,xl0,b,c,delt) 

implicit  real*8  (a-h,o-z) 

x  1 =f  1  case(x  1 0,-c,delt) 

f2cased  1  =(x20/  x  1 0+b/c*dlog(x  1 0/x  1 ))  *x  1 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccc 

real*8  function  f2cased3(x20,xl0,b,c,f,ddt) 

implicit  real*8  (a-h,o-z) 

xl=xlO+f*delt 

f2cased3=b/(f+c)*xl-(b*xl0/(f+c)-x20)*((xl0/xl)**(c/f)) 

return 

end 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccc 

real *8  function  fventfcn(pamb) 
implicit  real*8  (a-h,o-z) 

common/com  vent/pambO  1 0,pamb  1 03 0,  ventO  1 0,  vent 1 03 0,  vent30 

if(pamb.  lt.pambO  1 0)then 
fventfcn= ventO  10 
else  if(pamb.lt.pambl030)then 
f ventfcn= vent  1030 
else 

fventfcn=vent30 
end  if 
return 
end 


c 

End  of  source  code  for  program  MODEL.FOR 
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Sample  input  file  MODEL.INP  for  program  MODEL.FOR 


31,  !  number  of  nodes  in  profile 
0.85d0,  !  oxygen  consumption  rate  in  liter/min/diver 

lO.dO, !  co2  cone,  in  ppm  in  venting  gas 


6300.d0 

!  initial 

co2 

cone,  in  ppm  in  chamber 

O.dO, 

O.dO, 

O.dO, 

4, 

0, 

4.d0, 

O.dO, 

O.dO, 

4, 

1, 

5.18d0, 

O.dO, 

O.dO, 

4, 

0, 

8. dO, 

O.dO, 

4.d0, 

4, 

0, 

9.d0, 

O.dO, 

4.d0, 

4, 

1, 

10.18d0, 

O.dO, 

4.d0, 

4, 

0, 

12.d0, 

O.dO, 

4.d0, 

4, 

0, 

14.d0, 

20.d0, 

4.d0, 

4, 

0, 

15.d0, 

30.d0, 

4.d0, 

4, 

0, 

16.d0, 

30.d0, 

4.d0, 

4, 

0, 

21.d0, 

30.d0, 

4.d0, 

4, 

0, 

25. dO, 

30.d0, 

O.dO, 

4, 

1, 

26.d0, 

30.d0, 

O.dO, 

4, 

0, 

28.d0, 

30.d0, 

O.dO, 

4, 

0, 

31.d0, 

30.d0, 

O.dO, 

4, 

1, 

32.d0, 

30.d0, 

O.dO, 

4, 

o, 

33.d0, 

30.d0, 

O.dO, 

4, 

0, 

36.d0, 

30.d0, 

O.dO, 

4, 

1, 

37.d0, 

30.d0, 

O.dO, 

4, 

0, 

42.d0, 

O.dO, 

O.dO, 

4, 

1, 

43.18d0, 

O.dO, 

O.dO, 

4, 

0, 

47.d0, 

O.dO, 

O.dO, 

4, 

o. 

48.18d0, 

O.dO, 

O.dO, 

4, 

0, 

50.d0, 

O.dO, 

4.d0, 

4, 

o. 

53.d0, 

30.d0, 

4.d0, 

4, 

0, 

54.d0, 

30.d0, 

4.d0, 

4, 

0, 

55.d0, 

30.d0, 

4.d0, 

4, 

o, 

57.d0, 

30.d0, 

4.d0, 

4, 

o, 

58.d0, 

30.d0, 

4.d0, 

4, 

1, 

59.d0, 

30.d0, 

4.d0, 

4, 

0, 

60.d0, 

30.d0, 

4.d0, 

4, 

0, 

read(l,*)  nnode  !  number  of  nodes  in  profile 
read(l,*)  a  !  oxygen  consumption  rate 
read  (1,*)  cin  !  concentration  co2  in  ppm  in  venting  gas 
read(  1 , 1  l)tmin(i),fsw(i),dwater(i),ndiver(i),ivent(i) 

End  of  sample  input  file  MODEL.INP 
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The  next  several  pages  include  an  actual  input  data  file  (TM071797.INP)  that  was  used  for  the 
program  MODEL.FOR,  the  data  sheets  from  which  the  input  file  was  derived,  the  ouput  file 
(TO071797.OUT),  and  a  plot  of  the  calculated  CO2  levels  along  with  the  values  measured  using 
the  CO2  analyzer.  This  particular  data  set  was  chosen  due  to  its  richness.  The  person  recording 
the  data  did  an  excellent  job,  especially  in  recording  the  duration  of  the  venting  periods.  Note 
that  there  are  several  uncertainties  in  the  data:  (1)  The  depth  of  water  in  the  trunk  is  an  estimate, 
(2)  the  total  volume  of  the  unflooded  trunk  is  estimated  and  does  not  take  into  account  the  space 
taken  up  by  personnel  or  equipment  that  is  being  transfered,  (all  of  which  lead  to  uncertainty  in 
the  gas  volume),  and  (3)  the  concentration  of  CO2  in  the  venting  gas  is  unknown.  For  this 
particular  input  file,  CO2  was  assumed  to  be  present  at  a  very  low  level  of  10  ppm. 
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Datafile  TO071797.INP 


66,  !  number  of  nodes  in  profile 
0.85d0,  !  oxygen  consumption  rate  in  liter/min/diver 

lO.dO, !  co2  cone,  in  ppm  in  venting  gas 
4650.d0,  !  initial  co2  cone,  in  ppm  in  chamber 

0.,0.,0., 2, 1,4650, 0.77 
0.77,0.,  0.,2,0„ 

6..  0..0.. 2. 1.4340. 0.77 
6.77,0.,0.,2,0„ 

14.. 0..3..2.1. 3960.0.33 
14.33, 0.,3.,2,0„ 

21. . 0..3.5.2. 1.6040. 0.68 
21.68,0.,3.5,2,0„ 

29..  0.. 5.5. 2. 1.9220.0.77 
29.77,0.,5 .5,2,0,, 

36..  0.. 5.5. 2. 1.7820.0.65 
36.65, 0.,5.5,2,0„ 

43..  0..5.5.2. 1.7060.0.82 
43.82,0.,5.5,2,0„ 

49..  0.. 5.5. 2. 1.6640.0.78 
49.78,0., 5.5,2, 0„ 

57..  0.. 5.5. 2. 1.7280.0.75 
57.75, 0.,5.5,2,0„ 

64.. 5..5.5. 2. 1.7660.0.8 

64.8. 5. . 5.5. 2. 0„ 

72..  5.. 5.5.2. 1.7680. 0.52 

72.52. 5..  5.5. 2. 0„ 

79.. 5.. 5.5.2. 1.8200.0.82 

79.82. 5.. 5.5.2.0„ 

85..  5.. 5.5. 2. 1.6140.0.52 

85.52. 5..  5.5.2. 0„ 

92..  5..5.5. 2.1. 6420. 0.55 

92.55.5. . 5 .5.2.0., 

99..  5.. 5.5. 2. 1.7000. 0.97 

99.97. 5..  5.5. 2. 0„ 

105.. 5.. 5.5.2.1.4580. 0.63 

105.63. 5.. 5 .5.2.0., 

111..  5.. 5.5. 2.1. 5240. 0.77 

11 1.77. 5..  5.5.2. 0„ 

117.. 5.. 5.5.2. 1.4800. 0.63 

117.63. 5..  5.5.2. 0„ 

123.. 5..5.5. 2. 1.4840.0.73 

123.73.5.. 5 .5.2.0., 

129.. 5..5.5.2.1. 6880. 0.7 

129.7.5..  5.5.2.0„ 

136.. 5.. 5.5.2.1.5900. 0.92 

136.92.5..  5.5. 2. 0„ 

142. . 5. .5 .5.2. 1.3800.0.45 

142.45. 5..  5.5. 2. 0„ 

150..  5..5.5.2.1. 5700.1.05 
151.05, 5., 5.5, 2, 0„ 

156.. 5..5.5.2.1. 5200.0.67 

156.67.5..  5.5.2.0„ 
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164..  5..5.5. 2.1.7000.0.43 

164.43. 5..  5.5. 2.0„ 

170..  6..5.5. 2. 1.7400. 0.65 

170.65.6..  5.5. 2. 0„ 

175..  6..5.5. 2.1. 4720. 0.63 

175.63.6.. 5.5.2.0„ 

182.. 6..5.5. 2.1. 5840.0.58 

182.58.6.. 5.5.2.0„ 

189..  6.. 5.5. 2. 1.5600.0.47 

189.47. 6.. 5.5. 2. 0„ 

194..  6.. 5.5. 2.1.6000.0.58 

194.58. 6.. 5.5.2.0„ 

204..  6.. 5.5. 2. 1.5840. 0.58 

204.58.6.. 5.5.2.0„ 

215..  30.. 5.5. 2.0. 11380, 

216.. 40..5.5.2.0.13240, 

219. . 50..5 .5.2.0. 15160 

221.. 60..5 .5.2.0. 17780 

End  of  file  TO071797.INP 


26 


Output  data  file  TM071797.DAT 


********** 

.00000 

.00465 

1.00000 

4650.00 

.77000 

.00286 

1.00000 

2858.91 

6.00000 

.00425 

1.00000 

4247.10 

6.77000 

.00263 

1.00000 

2625.22 

14.00000 

.00496 

1.00000 

4964.50 

14.33000 

.00361 

1.00000 

3606.55 

21.00000 

.00649 

1.00000 

6485.02 

21.68000 

.00300 

1.00000 

3004.06 

29.00000 

.00852 

1.00000 

8516.59 

29.77000 

.00091 

1.00000 

914.85 

36.00000 

.00921 

1.00000 

9207.00 

36.65000 

.00126 

1.00000 

1264.69 

43.00000 

.00972 

1.00000 

9716.56 

43.82000 

.00089 

1.00000 

894.23 

49.00000 

.00779 

1.00000 

7788.83 

49.78000 

.00085 

1.00000 

850.64 

57.00000 

.01046 

1.00000 

10460.49 

57.75000 

.00109 

1.00000 

1089.71 

64.00000 

.00941 

1.15152 

9409.99 

64.80000 

.00121 

1.15152 

1206.31 

72.00000 

.01079 

1.15152 

10789.54 

72.52000 

.00253 

1.15152 

2528.59 

79.00000 

.01115 

1.15152 

11153.49 

79.82000 

.00130 

1.15152 

1300.19 

85.00000 

.00819 

1.15152 

8194.79 

85.52000 

.00201 

1.15152 

2005.67 

92.00000 

.01063 

1.15152 

10630.58 

92.55000 

.00232 

1.15152 

2315.33 

99.00000 

.01090 

1.15152 

10900.31 

99.97000 

.00097 

1.15152 

970.48 

105.00000 

.00767 

1.15152 

7665.42 

105.63000 

.00148 

1.15152 

1480.70 

111.00000 

.00863 

1.15152 

8628.19 

111.77000 

.00121 

1.15152 

1207.22 

117.00000 

.00817 

1.15152 

8168.36 

117.63000 

.00155 

1.15152 

1552.92 

123.00000 

.00870 

1.15152 

8700.41 

123.73000 

.00131 

1.15152 

1314.97 

129.00000 

.00833 

1.15152 

8329.36 

129.70000 

.00136 

1.15152 

1356.39 

136.00000 

.00974 

1.15152 

9741.71 

136.92000 

.00099 

1.15152 

990.11 

142.00000 

.00775 

1.15152 

7751.61 

142.45000 

.00227 

1.15152 

2270.77 

150.00000 

.01232 

1.15152 

12319.85 

151.05000 

.00091 

1.15152 

911.30 

156.00000 

.00750 

1.15152 

7499.77 

156.67000 

.00134 

1.15152 

1339.44 

164.00000 

.OHIO 

1.15152 

11095.69 

164.43000 

.00328 

1.15152 

3276.12 

170.00000 

.01069 

1.18182 

10690.11 

170.65000 

.00191 

1.18182 

1910.04 

175.00000 

.00770 

1.18182 

175.63000 

.00155 

1.18182 

182.00000 

.01003 

1.18182 

182.58000 

.00213 

1.18182 

189.00000 

.01068 

1.18182 

189.47000 

.00295 

1.18182 

194.00000 

.00898 

1.18182 

194.58000 

.00195 

1.18182 

204.00000 

.01449 

1.18182 

204.58000 

.00292 

1.18182 

215.00000 

.01679 

1.90909 

216.00000 

.01813 

2.21212 

219.00000 

.02212 

2.51515 

221.00000 

.02479 

2.81818 

********** 

End  of  output  file  TM071797.DAT 


7699.91 
1548.73 
10027.22 

2133.92 
10678.97 

2949.61 
8979.05 
1950.11 
14488.16 
2916.25 
16792.58 

18126.61 
22122.65 
24787.69 
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Following  are  the  two  raw  data  sheets,  as  well  as  a  graph  of  the  measured  and  calculated  C02  levels.  Bitmap  copies 
of  the  data  sheets  are  contained  in  files  DDSDATA1.PDF.  The  graph  is  contained  in  a  Hewlett-Packard  Graphics 
Language  file  TM071797.PDF. 
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